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We examine transport properties of a magnetic superlat- 
tice with current perpendicular to the planes. In the limit 
that the phase-breaking and spin flip scattering lengths are 
greater than the system size, a multiple-scattering approach 
is used to calculate the 4-probe conductance. We show that 
by tuning the strength of tunnel barriers placed between the 
current and voltage probes giant magnetoresistance ratios of 
arbitrary strength and size are achievable. 

Pacs numbers: 72.10.Fk, 72.15.Gd 



During the past few years advances in fabrication tech- 
nology have enabled the production of precisely engi- 
neered mesoscopic superlattice (SL) structures, consist- 
ing of alternating magnetic and non-magnetic metallic 
layers, with well defined dimensions and interfaces. Ex- 
change coupling of the magnetic layers through the non- 
magnetic material gives rise to antiferromagnetic (AF) 
alignment of adjacent magnetic layers. When such AF 
alignment is broken by applying a large magnetic field, a 
global ferromagnetic (F) configuration of the multilayer 
is achieved, and the resistance drops drastically 
Such functional magnetic materials are the focus of sub- 
stantial research both from a technological and funda- 
mental viewpoint, having a wide range of applications 
in magneto-electronics ||, such as read/write heads for 
high density magnetic storage systems, and in the minia- 
turization of magnetic field sensors, such as solid state 
compasses ||. All these applications require materi- 
als/geometries capable of giving a large signal (i.e. a large 
drop in resistance), and high sensitivity (so that only a 
small magnetic field required). Conventional materials, 
such as the 3d transition metal multilayers, show large 
GMR, but high magnetic fields are necessary to overcome 
the exchange coupling, making such systems unsuitable 
for most applications. For inhomogeneous multilayers 
based on polycrystalline alloys || much smaller fields are 
required, but show quite small drops in resistance. The 
half- metals such as Cr02, in principle can guarantee an 
infinite GMR @J|, but to date a GMR of only 50% has 
been measured in &O2 powders at low temperature || . 

In this letter we argue that attention should also focus 
on the problem of optimising the scattering properties 
of a given structure and on the role of external probes. 
By tunning the scattering potential in a 4 - probe geome- 
try we shall demonstrate a significant GMR enhancement 
can be achieved. 



Transport properties of GMR structures are typically 
studied, when the current flows either in the plane of 
the layers (CIP) or perpendicular to the plane of the 
layers (CPP). In the CIP configuration the dimensions 
are macroscopic and 4 probe measurement techniques are 
used [[l0| . In contrast CPP measurements are usually car- 
ried out with 2 probes, using sophisticated techniques to 
measure the small resistances involved (~ 10 -8 f2). The 
CPP configuration is particularly interesting because it 
produces a larger GMR signal, with changes in resis- 
tance up to 150% ]ri]Jl^| . Remarkably however very 
little research has been conducted using 4 probe mea- 
surements on CPP geometries |ll|, despite the fact that 
the quantum nature of transport in such structures gives 
rise to unexpected phenomena, as detected in anisotropic 
magnetoresistance (AMR) systems of Co/Ni multilayers 

00. 

In this letter we consider a 4 probe CPP structure, 
and demonstrate that due to an inherent instability of 
4 probe conductance measurements, an infinite magne- 
toresistance (IMR) is achievable. We use a general scat- 
tering approach to dc transport, developed to describe 
phase-coherent transport in dirty mesoscopic structures, 
and based on the fundamental current-voltage relations 
derived in [Lfj]. For simplicity in this Letter, we focus 
on the conductance of the structure shown in Fig. 1(a), 
which comprises a superlattice, with alternating normal 
and magnetic metal layers, in contact with four normal 
reservoirs at voltages Vj(j — 1 , . . . , 4) . The leads con- 
necting the structure to reservoirs 1 and 2 carry a cur- 
rent J, whereas leads 3 and 4 carry no current and hence 
form the voltage probes. The current/ voltage leads are 
separated by an insulating barrier. 

In what follows, we consider the zero-temperature, 
zero-bias limit in which the phase-breaking and spin-flip 
scattering lengths are greater than the dimensions of the 
system. In this limit, transport properties depend not 
only on the electronic structure of the magnetic multi- 
layers, but also on the contacts with external reservoirs. 
This feature is embodied in a fundamental current - volt- 
age relation due to Buttiker p6[, 
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which relates the current It from a normal reservoir i 
to the reservoir voltages (vj). The coefficients Gij sat- 
isfy Eli G -3 = Eli ^ = °- In units of 2e 7^ B 
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Gu — Ni — Ri and Gij-a — —Tij, where is the trans- 
mission coefficient from probe j to probe i, Ri is the 
reflection coefficient in probe i and Ni is the number 
of open scattering channels in lead i. As indicated in 
Fig. lb, Eq.(|l]) can be viewed as an equivalent circuit 
representation of the phase coherent structure of Fig. la. 
However the elements Gy are correlated functionals of 
the scattering potential generated by contact with exter- 
nal leads, geometry and disorder, and therefore cannot 
be varied independently. Furthermore, as noted in p^ ] 
Gij(H) = Gji(H*), where H is the Hamiltonian of the 
system and therefore in the presence of a magnetic field, 
Gij ^ Gji. 

Setting ii = — J 2 = /, h = h = and solving Eq.(||) 
for the conductance yields, 
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where d > is the determinant of the 3x3 matrix ob- 
tained by removing the third row and column from the 
Gij matrix of Eq.ffl). In what follows we compute the 
GMR ratio (G) as § 
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where Gj (Gf) is the conductance for the minority 
(majority) spin in the F configuration, and Gjj is the 
conductance for both spins in the AF configuration. 

In the F and AF configurations, the various transmis- 
sion and reflection coefficients for the separate spins can 
be computed by solving the Schrodinger equation on a 
nearest neighbor tight-binding lattice of sites. In sim- 
ple quasi one dimensional gemoetries involving only two 
probes, by combining an efficient recursive Green's func- 
tion approach with a material-specific s-p-d tight-binding 
Hamiltonian one can accurately predict the GMR ratio 
for a range of materials and layer thickneses pi] ]. For 
more complex geometries involving several probes, this 
is beyond the capabilities of currently available comput- 
ing resources and therefore to demonstrate a generic en- 
hancement of the GMR ratio, we analyze a tight binding 
model with a single degree of freedom per spin on each 
lattice site, as introduced in |]l7| , |l8| . Each lattice site is 
labeled by an index i and possesses a spin degree of free- 
dom tpf , and the corresponding Schrodinger equation has 
the form 



Eft = efVf - £>«&° 
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where j sums over all neighbors of i. The nearest 
neighbor hopping elements 7y fix the band-width and 
ef determines the band-filling. The separate spin fluids 
are assigned a spin-dependent on-site energy ef . The pa- 
rameters used in the present calculations are as follows: 



in the leads and non-magnetic regions of the SL ef — 1 
and 7y = 1; in the magnetic layers 7^ = 2 and ef = 1.45 
(3.7) for spin up (down). Following plj the hopping el- 
ements 7y joining magnetic-normal sites are chosen to 
be the geometric mean of the hopping elements in the 
magnetic/normal regions. 

The dc conductance (||) is determined by the multi- 
channel quantum mechanical scattering matrix, which 
in turn can be found from the Green's function for a 
given structure. The numerical approach involves at- 
taching the scattering region to crystalline semi-infinite 
leads, in which the Green's functions are superpositions 
of plane waves, and then using Dyson's equation to de- 
termine the total Green's function EJ. Effectively this 
procedure requires the inversion of Hamiltonian describ- 
ing the scatterer, which for the structure analyzed here is 
a 6080A6080 matrix per spin degree of freedom. As we 
are only interested in the Green's functions on the sur- 
face of the structure the inversion is carried out after the 
Hamiltonian has been renormalised in real space using 
an exact decimation technique [^3| . 

As shown in Fig. la the structure analyzed consists of 
four, crystalline semi-infinite leads each 40 sites wide. 
The superlattice consists of 5 non-magnetic layers alter- 
nating with 4 magnetic layers, each 8 sites wide and 40 
sites long. The leads are separated by barriers, 2 sites 
wide and of strength U. These barriers change the on- 
site energy ef in Eq.(Q) to ef + U. Hence, varying the 
magnitude of U alters the transmission and reflection co- 
efficients of the current/ voltage leads. For a barrier of 
width 2 sites, fig. 2 shows the dependence on U of the 
transmission coefficient per open channel, which as ex- 
pected, decays exponentially with U for large U . 

Fig. 3 shows the dependence of Gf,Gf j, and G| on bar- 
rier strength. We see that Gfj and G| increase smoothly 
as U is increased, whereas Gj diverges for U w 2.2, and 
then returns to a negative constant. This divergence oc- 
curs as the transmission coefficient from lead 4 to lead 
1 approaches the transmission coefficient from lead 4 to 
lead 2. 

To gain physical insight into the origin of this diver- 
gence, we return to the equivalent circuit of fig. lb and 
consider the case of a symmetric structure, where G13 = 
G24, G12 = G 34 , G14 = G23 and Eq.([| reduces to 

2(G43 + G42XG43 + G41) 
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G42 — G41 

From this expression we see that if G4 2 < G41 then 
the 4 probe conductance will be negative and when 
G42 = G41, the conductance is infinite. Of course since 
the quantities Gy are determined by the quantum me- 
chanical scattering properties of the structure in figure 
1, there is no guarantee that such a condition is achiev- 
able. The fact that we are able to demonstrate that this 
can occur in a realisable structure is a key result of this 
Letter. 
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Fig. 4 shows the dependence of the GMR ratio on bar- 
rier strength, calculated using Eq.(||). Due to the be- 
havior of the majority spin state (Gj), the GMR ratio 
increases with increasing U , diverges for U w 2.2 giving 
an IMR. Further increasing U causes the GMR ratio to 
return from -oo to a constant negative value, of around 
— 110%. We note that such negative GMR is not due 
to an intrinsic scattering asymmetry of the materials as 
measured in Ref. pj[ , but is due to the peculiar spin- 
dependent behavior of the 4 probe measurement. It is 
important to point out that the system we described is 
in principle tunable, for example by applying a gate elec- 
trode giving rise to the barrier potential U. Hence by 
changing the external gate potential, arbitrary GMR ra- 
tios may be obtained, using exactly the same principle 
as a Wheatstone bridge, despite the fact that the system 
described is fully quantum mechanical. 

In conclusion, we have shown for the first time that due 
to an inherent instability in the 4-probe conductance of 
mesoscopic structures an infinite GMR ratio is realizable. 
To-date a great deal of effort has been aimed at opti- 
mizing the materials used in GMR devices. The above 
results suggest that geometry and external gating may 
be equally crucial in the race to obtain inexpensive mag- 
netoresistive devices. This effect is a generic feature of 
the 4 probe geometry, and the use of a simple single-band 
tight-binding model is sufficient to demonstrate the possi- 
bility of IMR. Nevertheless more sophisticated material- 
specific tight-binding models [ pp|j2l| , or ab-initio calcula- 
tions Q , will be needed to provide quantitative insights 
into real device characteristics. 
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FIG. 1. Figure 1(a), The hashed area represents the mag- 
netic layers, and the clear areas represent normal regions. The 
black strip between the leads represent the portion of the vari- 
able-height barrier. Figure 1(b) shows the equivalent circuit 
of the structure shown in 1(a). 
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FIG. 2. The average transmission coefficient per open FIG. 4. The dependence of the GMR ratio on barrier 

channel for a 2D tight-binding lattice. The on-site energy is strength, 
eo and the barrier strength U. The barrier is 2 lattice sites 
wide, and the on-site energy is eo + U. 
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FIG. 3. The dependence of the 4 probe conductances 
GfjGfj. and Gi on barrier strength, for the structure shown 
in Fig.l. 
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